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Abstract. Using tools previously described and applied to the prototype galaxy NGC 891, we model the optical 
to far-infrared spectral energy distributions (SED) of four additional edge-on spiral galaxies, namely NGC 5907, 
NGC 4013, UGC 1082 and UGC 2048. Comparing the model predictions with IRAS and, where available, sub- 
millimeter and millimeter observations, we determine the respective roles of the old and young stellar populations 
in grain heating. In all cases, the young population dominates, with the contribution of the old stellar population 
being at most 40%, as previously found for NGC 891. After normalization to the disk area, the massive star- 
formation rate (SFR) derived using our SED modeling technique, which is primarily sensitive to the non-ionizing 
ultraviolet output from the young stellar population, lies in the range 7 x 10~* — 2 x 10~^ MQyr~^kpc~^. This 
is consistent with normalized SFRs derived for face-on galaxies of comparable surface gas densities from Ha 
observations. Though the most active star-forming galaxy of the five in absolute terms, NGC 891 is not an 
exceptional system in terms of its surface density in SFR. 
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1. Introduction 

Optical and Far-Infrared (FIR) /sub- millimeter (submm) 
data from galaxies contain complementary information 
about the distribution of stars and dust, from which intrin- 
sic quantities of interest - the star-formation rate (SFR) 
and star-formation history - can in principle be extracted. 
This is especially relevant to systems having intermediate 
optical depths to starlight, or to inhomogeneous systems 
with optically thin and thick components. Many, perhaps 
most, of normal (non-starburst) gas-rich galaxies in the 
local universe may fall into these categories. 

On the one hand, optical data probes the colour and 
spatial distribution (after correction for extinction) of the 
photospheric emission along sufhciently transparent lines 
of sight. This is particularly useful to investigate older, 
redder stellar populations in galaxian disks with scale 



heights larger than that of the dust. A radiation trans- 
fer modeling technique that can be applied to edge-on 
systems, where the scale height of the stars and dust ex- 
tinction can be directly constrained, was introduced by 
Kylafis & Bahcall ( IQSTj ) and subsequently applied on sev- 
eral galaxies (Xilouris et al. 1997: 1998; 199£, hereafter re- 
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ferred to as X97, X98, X99 respectively). Radiative trans- 
fer codes in combination with observations of nearby edge- 
on galaxies were also u sed by Ohta & Kodaira ( 1995 ) and 
Kuchinski et al. ( |l998D . 

On the other hand, grains act as test particles probing 
the strength and colour of ultraviolet (UV) / optical inter- 
stellar radiation fields. This constitutes an entirely com- 
plementary constraint to studies of photospheric emission. 
In the FIR, grains are moreover detectable over the full 
range of optical depths present in a galaxy. At least part of 
this regime is inaccessible to direct probes of starlight, es- 
pecially at shorter wavelengths, even for face-on systems. 
This particularly applies to light from young stars located 
in, or close by, the dust clouds from which they formed, 
since a certain fraction of the light is locally absorbed. 
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Furthermore, there is at least a possibility that most of the 
remaining UV and even blue light from young stars that 
can escape into the disk might be absorbed by diffuse dust 
there. Observations from IRAS, IRAM (e.g. Neininger et 
al. |1996|; Dumke et al. |1997D, SCUBA (e.g. Alton et al. 



1998| ; Israel et al. |1999| ; Dunne ^000| Bianc hi et a l. ^OOOaD 



and ISO (e.g. Haas et al. |1998| Stickel et al. |2000D can pro- 
vide information on the quantity and spatial distribution 
of the dust within spiral galaxies. 

A combined analysis of the whole UV-optical/FIR- 
/submm/mm spectral energy distribution (SED) of galax- 
ies seems to be a promising way to constrain the problem. 
Radiative transfer codes for an assumed "sandwich" con- 
figuration of dust and st ars w ere applied by Xu & Buat 
( |1995D and Xu & Helou ( |1996[ ) to account for the energy 
budget over the whole spectral range. They considered 
in detail the relative contribution of the non-ionizing UV 
photons and the optical photons in heating the grains. 
However these calculations did not incorporate a model 
for the dust grain emission, nor for radial variations in 
the absorbed radiation in the disk, and therefore could 
not account for the exact shape of the FIR SED. Recently, 
there have been several works modeling the SED of galax- 
ies from UV to submm (Silva et al. 199§| ; Devriendt et 
al. 1999), by applying photometric and/or spectrophoto- 
metric and chemical evolution models of galaxies. While 
these models are adequate in describing the general shape 
of the volume- integrated SED, they make use of many free 
parameters and of a simplified geometry. 

Bianchi et al. ( pOOObj ) attempted to model NGC 6946 
from the UV to FIR using a 3D Monte Carlo radiative- 
transfer code for a simplified geometry of emitters (a single 
stellar disk) . They concluded that the total FIR output is 
consistent with an optically thick solution. However, their 
model did not consider a grain size distribution for grains, 
stochastic emission of small grains, and the contribution 
of localized sources within star- forming complexes. This 
resulted in a poor fit of the FIR SED and a failure to 
reproduce the IRAS flux densities. 

Another work on modeling the UV to submm emis- 
sion was presented by Efstathiou et al. ( 200C| ) for star 
burst galaxies which were treated as an ensemble of op- 
tically thick giant molecular clouds centrally illuminated 
by recently formed stars. This modeling technique success- 
fully reproduced the observed SED of M82 and NGC 6090. 
Such a technique obviously cannot be applied to normal 
"quiescent" disk galaxies dominated by emission from the 
diffuse interstellar radiation field. 

In Popescu et al. ( 200C| , hereafter referred to as Paper 
I) we gave a detailed description of a new tool for a com- 
bined analysis of the optical-FIR/submm SED of galaxies. 
This tool included solving the radiative-transfer problem 
for a realistic distribution of absorbers and emitters and 
by considering realistic models for dust, taking into ac- 
count the grain size distribution and stochastic heating of 
small grains as well as the contribution of HII regions. We 
applied this tool to the edge-on system NGC 891. 



In short, we used the intrinsic spatial distribution of 
dust and stars (in the B, V, I, J, K bands) derived for this 
galaxy by X99 as a constraint on the old stellar popula- 
tion and part of the dust distribution. In order to fully 
explain the optical- FIR-submm SED, it was found neces- 
sary to add to the components of the X99's solution both a 
young stellar population (to correct for a shortfall in FIR 
brightness in terms of re-radiated non-ionizing UV) and 
more dust mass (to correct for a shortfall in submm bright- 
ness). Our solution explained both the observed SED in 
the FIR and submm, as well as the observed radial profile 
at 850 ^m. We found that the dust is predominanly heated 
by the young stellar population. 

Although NGC 891 was chosen for study in Paper I 
as the prototypical and most extensively observed edge-on 
galaxy, it exhibits some extreme features, raising the possi- 
bility that it might not in fact be representative of a "typ- 
ical" spiral. In this paper we extend our SED modeling 
technique to four additional edge-on systems - NGC 5907, 
NGC 4013, UGC 1082 and UGC 2048 - with the aim 
of examining whether the features of the solution we ob- 
tained for NGC 891 might be more generally applicable. 
Radiation transfer solutions accounting for the appear- 
ance of these four galaxies in the optical bands have been 
obtained by X99, X98 and X97. We use these solutions to 
constrain our model for the optical-FIR SEDs in the same 
fashion as for NGC 891. 

The paper is arranged as follows: In Sect. 2 we overview 
the model and its application to the objects in the present 
work. In Sect. 3 we present predictions for the FIR SED 
of all five galaxies studied. These are obtained using our 
"standard model", i.e. a model with radiation fields de- 
rived from the radiation transfer solution for the opti- 
cal appearance supplemented by a hidden population of 
young stars embedded in HII regions, but with no addi- 
tional dust. In Sect. 4 we describe and discuss the SED 
of one of our four galaxies, namely NGC 5907, in terms 
of the more complex "two-dust-disk model" and compare 
this with the corresponding solution obtained for NGC 891 
in Paper I. In Sect. 5 we discuss the dependence of the 
derived star-formation rates on the assumptions of the 
model, in particular the assumed geometries for the dust 
and stars. We show that after normalization to disk area, 
NGC 891 has comparable star-formation characteristics 
to the other four objects. We also compare the disk-area 
normalized massive star-formation rate in the five edge- 
on galaxies, obtained from the optical-FIR SED analysis, 
with the same quantity derived from H^ measurements of 
statistical samples of face-on galaxies given in the litera- 
ture. In Sect. 6 we give a summary of our work. 



2. Model calculations for the five galaxies 

A full description of our model and its application to 
NGC 891 can be found in Paper I. Here we give an 
overview of the basic idea and its specific application to 
the four galaxies. 
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To determine the propagation of light in a galaxy we 
take both emitters and absorbers to be smoothly arranged 
as a superposition of simple cylindrically symmetric geo- 
metrical distributions. The assumption of an axisymmet- 
ric model and the neglect of the spiral structure does not 
seem to introduce a systematic error on the estimate of 
the overall opacity, at least when a statistical sample of 



edge-on galaxies is considered (Misiriotis et al. 200C). The 
emitters are divided into a young (predominantly UV- 
emitting) stellar population, distributed in a disk with 
a small scale height, an old (predominantly optical-NIR 
emitting) population, distributed in a disk with a larger 
scale height and a bulge. The absorbers are arranged in 
"young" and "old dust disks" , associated with the old 
and young stellar populations. The "old dust disk" has 
a larger scale height than the "young dust disk" , though 
not so large as the old stellar disk. Both the stellar- and 
dust-disk distributions are specified as exponential density 
distributions in radius and height, while the stellar bulge 
is described by a de Vaucoulcurs law. All radial distribu- 
tions have a common truncation at three scale lengths of 
the "old dust disk" . 

The first step in the calculation of the optical-FIR SED 
is to determine the geometric and amplitude parameters 
for the old stellar disk and bulge independently in each 
band pass ranging from the B band to the NIR for a com- 
mon "old dust disk" (whose geometry and opacity is also 
determined). This can be achieved by a radiation trans- 
fer analysis matching as closely as possible the predicted 
and observed brightness distributions, following Kylafis & 



Bahcall (1987) and including inclination as a free param- 
eter. For the five galaxies discussed in this work we adopt 
the solutions given in X97 and X99 on the basis of obser- 
vations in B, V and I, and in the case of NGC 891 also in 
J and K bands. 

This analysis of the optical-NIR emission also yields 
the extinction coefficient as a function of wavelength for 
the "old dust disk" . For all five galaxies studied, the mea- 
sured extinction law was found to be close to the predicted 
law derived from the dust model used in the calculation of 
the IR emission. This corresponds to the graphite-silicate 
mix and dust size distribution determined for the Galaxy 
by Mathis et al. ( |1977D (see Eqs. [2] and [3] of P aper I ) 
with optical constants taken from Laor & Draine ( 1993| ). 
When calculating the IR emission we take grains in the 
"young dust disk", if present, to have the same size dis- 
tribution and optical properties as those in the "old dust 
disk". 

To this fixed basis for the old stellar population and 
associated dust disk, the model allows the addition of a 
young, UV-emitting stellar population and, optionally, as- 
sociated dust. The additional stars and dust are taken to 
be distributed in a common exponential disk with a scale- 
length equal to that of the already determined intrinsic 
B-band population, and a scale height fixed at 90 pc (as 
for the Milky Way). Since this scale height is typically 
a factor of several times smaller than the scale heights 
derived by X97 and X99 for the "old dust disk" , this con- 



stitutes the simplest assumption for the distribution of 
these extra components, which hide UV-optical indicators 
of the young stellar population and associated additional 
opacity. As emphasized in Paper I, a clumpy distribution 
for the dust associated with the young stellar population 
is also possible, though very difficult to calculate, and, due 
to the lack of high angular resolution FIR images of spiral 
galaxies, at present not directly verifiable. The luminosity 
of the young stars is a free parameter of the model, which 
we express as a recent star-formation rate SFR, based 
on the population synthesis models of Bruzual & Chariot 
( 2001 ) for Z — Zq, a Salpeter IMF with a cut-off mass of 
100 Mq and SFR = (l/r)exp(-i/T) with r = 5Gyr. 

The model also incorporates a parameter F represent- 
ing the fraction of the non-ionizing UV luminosity from 
the young stellar population, which is locally absorbed 
in star-formation regions by optically thick fragments of 
their parent molecular clouds. This fraction of the lumi- 
nosity is considered to be re-radiated in the FIR according 
to a spectral template matching the observed spectrum of 
galactic HII regions. This template generally has a warmer 
60/100 /im colour than that calculated for the diffuse disk. 

Given the intrinsic distributions of stellar emissivities 
in the UV, optical and NIR, and the distribution of ab- 
sorbers in the "old" and (if present) the "young dust 
disks" , we then proceed with a radiative transfer calcu- 
lation to determine the UV-NIR energy density of the 
radiation field throughout the galaxy. In the absence of 
solutions in the J and K bands, we extrapolate the energy 
density from the B, V, I bands assuming a black-body 
spectral distribution and a colour temperature of 4000 K. 
In calculating the heating of grains placed in the resulting 
radiation field, we include an explicit treatment of stochas- 
tic emission as detailed in Paper I. The IR-submm emis- 
sion from grains is then calculated for a grid of positions 
in the galaxy. Subsequently we integrate over the entire 
galaxy to obtain the IR-submm SED of the diffuse disk 
emission. Prior to comparison with observed FIR-submm 
SEDs, the spectral template for the HII regions, scaled ac- 
cording to the value of F, must be added to this calculated 
spectral distribution of diffuse FIR emission. 

Due to the precise constraints on the distribution of 
stellar emissivity in the optical-near infrared (NIR) and 
the distribution of dust in the "old dust disk" yielded 
by the radiation transfer analysis of the highly resolved 
optical-NIR images, coupled with the simple assumptions 
for the distribution of the young stellar population and 
associated dust, our model has a maximum of three free 
parameters from which the FIR-submm SED can be fully 
determined. These are the SFR for massive stars, the pa- 
rameter F representing the fraction of the non-ionizing 
UV luminosity from massive stars which is locally ab- 
sorbed, and the mass of dust in the "young dust disk". 
In general terms, the last parameter is principally con- 
strained by the submm emission, the non-ionizing UV lu- 
minosity by the bolometric FIR-submm output and the 
factor F (in the absence of high resolution images) by the 
FIR colour. In the analysis of the five galaxies in subse- 
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quent sections, we refer to the following variants of the 
model: 

1. The simple model: The simplest possible model takes 
into account only the information from the optical ob- 
servations, i.e. the dust and old stellar distribution in 
the old disk without massive UV-emitting stars. The 
sole purpose of the "simple model" is to demonstrate 
that the old stellar population can not adequately heat 
the dust to account for the observed FIR luminosity. 

2. The standard model: This is the "simple model" sup- 
plemented by a disk population of UV-emitting stars 
without associated dust, but with localised FIR emis- 
sion. The free parameters are the SFR and the factor 
F. 

3. The two-dust-disk model: This is the "standard model" 
supplemented by the presence of dust in the young stel- 
lar disk. The free parameters are the SFR, the factor 
F, and the dust mass in the "young dust disk" . 

The predictions of the model calculations are com- 
pared with available FIR, submm and mm data from the 
literature. In Table 1 we summarise the data used for com- 
parison in the present work. We adopted 20% errors in the 
IRAS and IRAM data. For the SC UBA d ata we adopted 
the errors quoted by Alton et al. (1998). In Table 1 we 



also give some basic information on the properties of the 
sample galaxies. One factor in the choice of these galaxies 
for the radiation transfer analysis by X97 and X99 was 
their prominent dust lanes. 

3. Results for the "standard model" 

To demonstrate the need for a population of massive stars 
to account for the observed FIR luminosities in all five 
galaxies, we first calculated SEDs for the "simple model" , 
in which the stellar and dust distributions are completely 
determined from the optical-NIR data as described in 
Sect. 2. These are plotted in the left column of Fig. |] 
In all cases, the predicted FIR flux densities in the 60 and 
100 /xm IRAS bands fall short of the observed luminosities 
by factors of between 5 and 10, showing that the heating 
of the dust from the old stellar population cannot account 
for the FIR emission. The results of the "simple model" 
are summarised in the second column of Table ^ (^simple), 
where we give the FIR luminosity of the galaxies that can 
be attributed to the old stellar population. 

Assuming that the optically determined dust content 
we used for each galaxy is exact, the "standard model" can 
account for the observed FIR luminosity by transforming 
non-ionizing UV radiation from massive stars into the FIR 
through a combination of local absorption and absorption 
in the diffuse old disk, controlled by the factor F. If all the 
young stellar luminosity is absorbed locally {F — 1), then 
the efficiency of the transformation from non-ionizing UV 
to the FIR is unity, as for starburst galaxies. In the total 
absence of local absorption {F = 0) the global absorption 
efficiency is lower, with a value depending on the relative 



geometries of the diffuse dust and the massive stars. In 
principle, the factor F can be determined from the FIR 
colours. 

We investigated the combinations of F and SFR for 
which the "standard model" was consistent with the 
IRAS 60 and 100 ^m measurements. For NGC 891 and 
UGC 2048, good agreement between model and data 
could be found for values of F well above zero. However, 
for the remaining galaxies, the IRAS colours required F 
values of zero or close to zero. This would imply high 
current star-formation rates and an amplitude in the 
UV stellar emissivity which would not join up smoothly 
with the empirically determined B band stellar emissivi- 
ties. Consequently, we have plotted in Fig. the model 
prediction corresponding to the maximum value of F 
still consistent with the IRAS colours. This also gives a 
monotonically increasing SED for the UV-optical stellar 
emissivities. Corresponding FIR bolometric luminosities 
(-^standard), SFRs and factors F are given in Table ||, 
columns 3, 4, 5 respectively. We emphasise however, that 
in principle one could have a non-steady SFR, so that the 
tabulated SFR are at the same time lower limits in SFR 
and upper limits for the factor F. Therefore we can dis- 
cuss the results for the "standard model" also in terms of 
a range of possible values for F and SFR, which will be 
compatible with the range in the IRAS errors. 

We remark that a comparison of the values in the sec- 
ond and third column of Table ^ reveals that, according to 
the standard model, only 10% to 30% of the dust emission 
is powered by the old stellar population. 

In the above we have made the assumption that the 
"standard model" is a sufficiently realistic model of the 
galaxies. Some evidence that this may not be the case is 
apparent from Fig. ^ (right column), where all the model 
spectra peak at around 120 fim, well shortward of the typ- 
ical peak near 150 /xm seen in spiral galaxies. This may be 
an indication that the overall optical depth of the galax- 
ies is underestimated by the "standard model" , indicating 
the presence of a further dust component in the plane 
of each galaxy. This component may either be smoothly 
distributed or be in the form of optically thick clumps 
or some combination of these. For the two galaxies that 
submm and/or mm measurements are available (NGC 891 
and NGC 5907) , this is explicitly confirmed by the under- 
shoot of the "standard model" prediction in the submm 
and mm regime. In Paper I, we used the available submm 
ffuxes to quantify the missing dust and introduced the 
"two-dust-disk model" , which successfully reproduced the 
spectrum and the spatial distribution of the FIR emission 
in NGC 891. In the next section we will apply the "two- 
dust-disk model" to NGC 5907 to determine the total dust 
content of this galaxy and examine the effect of this more 
realistic case on the derived SFRs. 
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Fig. 1. The SED of the 5 galaxies. In the left column the SED is derived for the "simple model". In the right column 
the SED is derived for the "standard model" . The predicted FIR emission from the diffuse disk is plotted with dashed 
lines and the contribution of localised sources with dotted lines. The corresponding SFR and factor F for each galaxy 
are given m Table |. For NGC 891 and UGC 2048, the model corresponding to the value of F which best fits the IRAS 
colour is plotted. In the cases of NGC 5907, NGC 4013 and UGC 1082, the model corresponding to the maximum 
value of F still consistent with the IRAS colour is plotted (see text). 
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Table 1. Properties of the galaxies used in the modeling of the FIR-submni SEDs. The dust mass is that quoted in 
X97 and X99 (corresponding to the "old dust disk" in this paper). The submm fluxes are taken from ^ Alton et al. 
( |1998D , ^Young et a l. ( |1989D , ^Soifer et al. ( |1989D , ^Moshir et al. ( |1990[ ). The mm fluxes from ^Dumke et al. ( [1997D , 
^Guelin et al. p9^ . 



Galaxy 


Hubble 


distance 


dust mass 




-F'loo 


-F450 


Fsoo 


-F1200 


-F1300 




type 


(Mpc) 


(Mq X lO'^) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


NGC 891 


Sb 


9.5 


5.6 


50.5" 


126' 


32' 


4.62' 




0.85'' 


NGC 5907 


Sc 


11.0 


1.5 


16.3^ 


55.9^ 






0.54^ 




NGC 4013 


Sbc 


11.6 


0.45 


7.0^ 


23.1^ 










UGC 1082 


Sb 


37.0 


0.99 


1.6* 


3.7* 










UGC 2048 


Sb 


31.5 


3.5 


1.7* 


3.5* 











' The flux densities from this reference are integrated in the range ±225" along the major axis of the galaxy. 



Table 2. Results for the "standard model". Lsimpic de- 
notes the bolometric re-radiated luminosity when only the 
old stellar population (B, V, I, J, K bands) is taken into 
account (i.e. for the "simple model"), ^standard denotes the 
bolometric re-radiated luminosity obtained by integrating 
the solution for the "standard model" as plotted in Fig. |l|. 
Columns 4 and 5 contain the corresponding SFRs and fac- 
tors F. 



Galaxy 


-^'simple 


-^/standard 


SFR 


F 




WxlO''^ 


WxlO^^ 


Mgyr-' 




NGC 891 


10.8 


48.1 


3.48 


0.28 


NGC 5907 


3.7 


15.5 


2.17 


0.08 


NGC 4013 


2.0 


11.0 


1.23 


0.19 


UGC 1082 


2.3 


21.3 


2.59 


0.23 


UGC 2048 


5.7 


18.3 


1.36 


0.27 



4. Application of the "two-dust-disk model" to 
NGC 5907 

In Paper I we proved that in NGC 891 the amount of 
dust (on the basis of the grain optical constants adopted) 
quoted by X99 is underestimated by a factor of two, and 
that the addition of a second dust disk, that follows the 
spatial distribution of the young stellar population, repro- 
duces very well the observed FIR luminosity and the radial 
profiles observed by Alton et al. (|l998D . For NGC 5 907, 
the measurement at 1200 /xm by Dumke et al. ( 1997 ) can 
be used to constrain the quantity of dust in the second 
disk. 

Inspection of Fig. |l| (right column) shows that the 
"standard model" solution for F = 0.08, which is in 
agreement with the IRAS points, crassly undershoots the 
1200 /im point. As for NGC 891, the only way to account 
for both the IRAS and the 1200 /im observations is to in- 
clude more dust in our model, as simply increasing the 
star-formation rate would provide too much luminosity in 
the IRAS range. 

The second-dust-disk's spatial distribution was con- 
strained to follow that of the young stellar population, 
which is taken to be exponential with scalelength and 
scaleheight = 5020 pc and — 90 pc respectively (see 
Sect. 2). To find the solution, the total amount of dust 
in the second dust disk was varied jointly with the SFR 



until the FIR-submm SED of our model fitted all three 
measurements. In Fig. |^ we present the resulting spec- 
trum overlaid with the data. This spectrum corresponds 
to a SFR= 2.2 MQyr^^ and F = 0.10. The total dust mass 
in the second dust disk is 4.5 x 10^ Mq. By comparison, 
the dust mass quoted for the old disk of NGC 5907 by 
X99 is 1.5 X 10"^ Mq. 

Thus, we have had to add 3 times more dust over 
the quantity implied by the observed dust lane to obtain 
agreement with the data on the basis of the "two-dust-disk 
model" . The derived value of F is still very low, but the 
spectrum of the diffuse component is too warm (NGC 5907 
is more optically thin than NGC 891) to allow the addi- 
tion of a warm component related to HII regions. While 
we could have produced a colder spectrum for the diffuse 
emission, allowing some room for HII regions, this would 
have been at the cost of adding yet more dust in the sec- 
ond disk. A larger quantity of dust in the second disk does 
not seem very likely, as the total dust mass for the fitted 
F = 0.10 case already implies a gas-to-dust ratio of 130 by 
mass (instead of 520 according to X99). As a comparison, 
the gas-to-dust ratio in the Milky Way is 133. 

The total IR-submm re-radiated luminosity of 
NGC 5907, obtained by integrating the F = 0.10 "two- 
dust-disk model" SED, is 50.5 x lO^^W out of which 
27.0 X 10'^^ W is attributed to heating from the young stel- 
lar population. Thus, about 40% of the dust emission is 
powered by the old stellar population. More dust in the 
"young dust disk", a lower SFR and a higher factor F 
would increase the disk opacity and thus the contribution 
of the old stellar population, but we excluded this solution 
due to the high dust-to-gas ratio. The 40% contribution 
from the old stellar population represents an upper limit, 
since a possible dumpiness in the dust distribution would 
decrease the absorption efficiency in the optical band. 

On the basis of the model's inherent assumption, that 
all the dust (in both disks) of the galaxy is distributed 
smoothly, the central face-on optical depth in the optical 
band is Tv = 1.4 for NGC 5907. 

The major difference between NGC 891 and 
NGC 5907, on the basis of the "two-dust-disk model" , 
is that the spectrum of the former apparently allows for 
the existence of a larger contribution from HII regions. 
Presumably NGC 5907 must have HII regions, and there 
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Fig. 2. SED of the "two-dust-disk model" for NGC 5907 
for F = 0.10 and SFR = 2.2M0yr-\ overlaid with data 
points given in Table 1. The predicted FIR emission from 
the diffuse disk is plotted with a dashed line and the con- 
tribution of localized sources with a dotted line. 

is no obvious reason why the local properties of these HII 
regions, in particular the probability of local absorption of 
the UV radiation, should differ between the two systems. 

One potential way out of this difficulty might be to 
suppose that the spectral template for the HII regions used 
in our model underestimates the associated submm emis- 
sion. As stated in Paper I, we did not attempt to include 
potential cold dust emission components that might be 
expected from "parent" molecular clouds in juxtaposition 
to their "offspring" HII regions. These might contribute a 
significant fraction of the observed submm emission, but 
be relatively faint in the IRAS range. If so, the dust con- 
tent and the luminosity of the diffuse disk would be lower 
for both galaxies, but more so for the more optically thin 
diffuse disk of NGC 5907. This would make room for larger 
values of F in the model fits, particularly for NGC 5907. 

An alternative scenario, also discussed in Paper I, 
would be the presence of truly quiescent cold dust clouds 
in the disk not associated with HII regions. In that case, 
the difference between the two galaxies might indicate a 
greater star- formation activity in the disk of NGC 891, in 
the sense that a larger proportion of the cold dust clouds 
in the disk had been triggered into massive star formation. 

5. Discussion 

Our results confirm that it is a common tendency for 
edge-on galaxies to hide a significant fraction of their dust 
from optical extinction studies. In two cases where sub- 
or near- mm data are available (NGC 891, as analyzed in 
Paper I and NGC 5907 in this work), we have modeled 
the FIR-submm SED assuming that the additional dust 



is distributed in a second diffuse disk associated with the 
young stellar population. We have found that the SEDs 
could be fitted in terms of total dust masses of respec- 
tively twice and four times the masses inferred by X99 
from the optical radiation transfer analysis for NGC 891 
and NGC 5907. The corresponding SFRs and local non- 
ionizing UV absorption factors F in the best fits were 
SFR = 3.8M0yr-i and F = 0.22 for NGC 891 and SFR 
= 2.2M0yr-i and F = 0.10 for NGC 5907. 

The main uncertainty in the inferred SFRs is the value 
of the factor F. In principle, F can be determined from the 
colour of the FIR SED, especially the 60/100 /^m colour 
ratio. However, as discussed in Sect. 4, the value of the fac- 
tor F, which can be fitted, in practice is constrained by 
the assumption of the "two-dust-disk model" that all the 
additional dust needed to fit the submm measurements 
is diffuse. In reality, the geometry of the additional dust 
could be, as emphasised in Paper I, in the form of optically 
thick clumps associated with the HII regions. In this case 
the inferred SFRs would, depending on the factor F, differ 
from those determined on the assumption of a purely dif- 
fuse dust distribution without local absorption. However, 
at present, the available data does not allow us to quantify 
the clumpy component of UV- absorbing dust. 

The SFRs and factors F inferred from the "standard 
model" for NGC 891 and NGC 5907 in Sect. 3 (where the 
dust content and distribution is determined from X99's 
optical analysis) hardly differs from the values determined 
from the "two-dust-disk model" applied to NGC 891 and 
NGC 5907 in Sect. 4. This is at first sight surprising, as 
one would expect lower SFRs due to the higher opacity 
of the diffuse disk in the "two-dust-disk model" . However, 
the SFR determined from the "standard model" was ob- 
tained by fitting the SEDs exclusively to the IRAS 60 
and 100 fim measurements common to the sample of the 
5 galaxies. The reason that the SFRs determined from 
the "standard model" are so similar to those derived from 
the "two-dust-disk model" is simply that the additional 
submm luminosity arises from the extra dust (from the 
second dust disk) rather than from additional SFR. As re- 
marked in Sect. 3, the rather warm SEDs for all the "stan- 
dard model" fits suggest that additional dust is present in 
the disks of all 5 galaxies. 

To statistically evaluate our results for SFR, we make 
the assumption that the SFRs determined from the "stan- 
dard model" fits to the IRAS 60 and 100 /im measurements 
for our sample of 5 galaxies give, as is the case of NGC 891 
and NCG 5907, a reasonable estimate of the SFR which 
would have been obtained with an analysis based on the 
"two-dust-disk model". We can then compare the SFR 
characteristics of the 5 edge-on galaxies with the larger 
sample of 61 galaxies with inclinations less than 75 de- 
grees which were studied by Kennicutt ( 1998] ) on the basis 
of Hq, measurements. 

Fig. ^ depicts the relation between the disk-averaged 
surface density in SFR (Ssfr) determined from the "stan- 
dard model" as a function of the average gas surface den- 
sity (Sg) for the 5 galaxies. The plotted data is summa- 
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rized in Table ^. The Ssfr is calculated from the SFRs de- 
termined from the "standard model" and listed in Table ^. 
The full range in conceivable SFRs is given by the verti- 
cal error bars. The upper and lower limits for the SFRs 
are calculated such that the predicted SED is still consis- 
tent with the IRAS colours (within the 20% IRAS error 
bars). Lower error bars are not given when the plotted 
SFRs represent lower limits for the SFR (maximum val- 
ues for the factors F, see discussion in Sect. 3). In these 
cases lower limits would be possible only if we allowed 
for different sources of uncertainties, like variations in the 
spectral shape of the template used for the HII regions. 
However this is hard to quantify, and in the following we 
assume that the errors of the SFR are given only by the 
uncertainties in the IRAS data. 

The Eg were calculated from the gas mass es taken from 
Garcia-Burillo et al. ( 1992 ) and Rupen ( 1991 ) respec tively 
for H2 and HI in NGC 891, from Dumke et al. ([1997D 
for H2 and HI in NGC 5907, from Bottema ( |l99q ) and 
Gomez & Garcia ( |1997D for H2 and HI in NGC 4013, from 
Giovanelli & Haynes ( |199^) fo r HI in UGC 1082 and from 
Huchtmeier & Richter ( |1989| ) for HI UGC 2048. In the 
last two cases, where only data for HI was available, we 
estimated the total gas mass by doubling the HI gas mass. 
The horizontal error bar corresponds to the uncertainty in 
the gas masses for which we adopted an average 0.2 dex 
error. The surface area of the disk was calculated for Ro = 
(3 ± 0.5)/id, where ha is the intrinsic radial scalelength 
determined from the radiation-transfer modeling in the I 
band. In their analysis of surface photometry of the outer 



regions of spiral disks, Pohlen et al. (2000) show that the 
disk boundaries are typically in this range. 

The points for the 5 galaxies in Fig. |^ lie within 
the area of the diagram occupied by the galaxies in the 
Kennicutt sample. The match is even better for those 
members of the Kennicutt sample with Hubble types Sb 
to Sc. This agreement is quite reassuring, bearing in mind 
the several factors which could introduce a systematic dif- 
ference between the SFRs inferred for a sample of nearly 
face-on systems from H^ measurements compared with 
the present technique for edge-on systems based on an 
analysis of broad-band non-ionizing UV re-radiated in the 
FIR-submm range. 

Firstly, the H^ analysis is sensitive to the most massive 
stars and in particular to the assumed mass cut of 100 M©. 
Whereas the FIR-submm modeling also assumes the same 
mass cut in the conversion of SFR to non-ionizing lumi- 
nosity (see Sect. 2 and Paper I), our model is less sensitive 
to this effect. 

Secondly, whereas the H^ is sensitive to the star- 
formation history of the last lO^yr, our broad-band FIR- 
submm SED analysis samples approximately the last 
lO^yr. Thus, our analysi s is c onsistent with the basic hy- 



pothesis (see Kennicutt 1998 ) for "normal" spiral galax- 
ies of a steady star-formation activity. In principle, we 
could extend our analysis based on our determinations of 
the intrinsic populations and use the determined intrinsic 



Table 3. The average gas surface density (Sg) and disk- 
averaged SFR surface density (Esfr) for our galaxy sam- 
ple. The lower and upper limits in the Esfr are also given 
together with the disk area used to normalize the SFRs 
and gas masses. 



Galaxy 



log Eg 

Mo 



log EsFR 

M0 



logEgiJ'A 

Mo 



yr kpc^ yr kpc^ 



logEgV^^^ 
Mo 
yr kpc^ 



Area 
kpc^ 



NGC 891 


1.09 


-2.30 


-2.54 


-2.07 


687 


NGC 5907 


1.27 


-2.29 


-2.29 


-2.15 


421 


NGC 4013 


1.40 


-1.85 


-1.85 


-1.63 


88 


UGC 1082 


0.87 


-2.27 


-2.27 


-1.89 


482 


UGC 2048 


0.70 


-2.80 


-3.17 


-2.60 


855 



colours to determine more accurately the SFR history of 
the galaxies, though this is beyond the scope of this paper. 

The assumption of a steady-state star-formation rate 
is also broadly consistent with the timescales for the ex- 
haustion of the current gas supply under the derived SFRs. 
The dotted, dashed and dot-dashed lines in Fig. |^ repre- 
sent star-formation efficiencies corresponding to consump- 
tions of 100, 10 and 1 percent of the gas mass in 10^ yr. 
We note that in constructing Fig. 3, we have assumed that 
the HI line emission, from which the gas masses for the 5 
edge-on galaxies were in part determined, have not been 
affected by self absorption in the 21cm line. The effect of 
self- absorption would be to shift the points for the 5 edge- 
on galaxies to the left of the positions corresponding to 
their actual gas surface densities in Fig. ^. 

Thirdly, the SFRs derived from Hq, were corrected by 
a single factor for extinction, despite the varying orien- 
tations. As well as possibly affecting the vertical position 
of the galaxies on the plot, this may induce some scatter, 
especially if all the dust were diffusely distributed. The 
systematic effect may be expressed in terms of the factor 
F: an overestimation of the factor F is equivalent to an 
overestimation of the local extinction in the star-formation 
regions (statistically averaged over the population of HII 
regions in a disk). Thus, while moving to higher factors 
F would move the points for the 5 galaxies towards lower 
SFRs, it would have the opposite effect for the SFRs de- 
termined from the Hq. 

Lastly, we remark that NGC 891 does not appear as 
an exceptional system compared with the other 4 galaxies 
in our sample (and with Kennicutt's [1998] normal galaxy 
sample) on the basis of SFR normalized to disk area. Our 
work thus provides no evidence that this galaxy's excep- 
tional layer of extraplanar He-emitting diffuse ionizing gas 



(e.g. Hoopes et al. 1999) and surrounding X-ray-emitting 
hot gas (e.g. Bregman & Houck 1997) is attributable to 
unusual star-formation activity. 

6. Summary 

— This paper is the second part on a series of papers 
dedicated to the modeling of the SEDs of disk galax- 
ies. In Paper I we described a new tool for the analysis 
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Fig. 3. Disk-averaged SFR surface density (Ssfr) as a 
function of average gas surface density (Sg) for our galaxy 
sample and for the sample of Kennicutt (1998) of 61 nor- 
mal disk galaxies with SFR determined from measure- 
ments. The 5 galaxies from our sample are plotted as filled 
circles and the corresponding Ssfr, 5]™pp, S™p^, Sg and 
disk areas are listed in Table |^. The SFR surface densities 
were calculated by averaging the SFRs determined from 
the "standard model" (Sect. 3) over disks with an optically 
defined boundary (i?o) taken to be 3 times the intrinsic 
radial scalelength hs determined from the radiation trans- 
fer modeling in the I band. The galaxies from the sample 
of Kennicutt are plotted as open circles (Sb,Sc,SBb,SBc), 
triangles (Sa), open squares (Unknown/Not Available), 
crosses (Irr). The dotted, dashed and dot-dashed lines 
represent star-formation efficiencies corresponding to con- 
sumptions of 100, 10 and 1 percent of the gas mass in 
10^ yr. 



of the UV to the submm SED and applied this tool 
to the well known nearby edge-on galaxy NGC 891. In 
the present paper we have extended the analysis from 
Paper I to a sample of four additional edge-on galaxies 
(NGC 5907, NGC 4013, UGC 1082, UGC 2048) and 
proved that the solution obtained for NGC 891 is gen- 
erally valid. We have also shown that NGC 891 is not 
an exceptional system in terms of its SFR density. 
In all the systems the dust is predominantly heated by 
the young stellar population. The contribution of the 
old stellar population can account for at most 40%. 
The SFRs derived from our "standard model" are 
consistent with the SFRs derived from the more so- 
phisticated "two-dust-disk model" and are also consis- 
tent with the Kennicutt 's (|1998D Schmidt law for disk 
galaxies. 
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